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The production of neutral pions in 158-A GeV ^"^Pb+^^^Pb collisions has been studied in the 
WA98 experiment at the CERN SPS. Transverse momentum spectra are studied for the range 
0.3GeV/c < mT — mo < 4.0GeV/c. The results for central collisions are compared to various 
models. The centrality dependence of the neutral pion spectral shape and yield is investigated. An 
invariance of the spectral shape and a simple scaling of the yield with the number of participating 
nucleons is observed for centralities with greater than about 30 participating nucleons which is most 
naturally explained by assuming an equilibrated system. 
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Ultra-relativistic heavy-ion collisions produce dense 
matter which is expected to be in the form of a de- 
confined phase of quarks and gluons at sufficiently high 
energy densities. The transverse momentum spectra of 
produced pions can provide information on both the ini- 
tial and final state properties of the hot hadronic matter. 
The high pT pion production is expected to be dominated 
by hard scattering of the partons. In pA collisions the 
high pt region is known to be enhanced (Cronin effect 
) due to initial state scattering of the incident partons 
leading to a broadening of their incoming px- In AA col- 
lisions, many of the scattered partons must traverse the 
excited matter to escape and therefore may undergo ad- 
ditional rescatterings and energy loss In the case of 
significant parton rescattering, the parton distributions 
may approach thermal distributions with a temperature 
reflecting the initial state of the excited matter. The in- 
termediate Pt region of the pion spectrum might then 
reflect this initial temperature while the low px region 
would dominantly reflect the freeze-out temperature. In- 
deed, one of the earliest signatures of QGP formation, 
proposed by Van Hove |^ , was the observation of a satu- 
ration of the average transverse momentum with increas- 
ing energy (or entropy) density for systems excited just 
above the critical energy density. With increasing energy 
density, the initial temperature would not rise above the 
critical temperature until all of the latent heat of the 
QQP phase transition had been extracted. 

For these reasons it is of interest to study the central- 
ity dependence of the pion production. It is generally 
believed that the initial energy density increases with in- 
creasing centrality, due to the many overlapping interac- 
tions. Also, the volume of the excited matter increases 
with centrality, as well as the amount of rescattering. 
Since rescattering is the feature which distinguishes AA 
collisions non-trivially from pp collisions, and since sig- 
nificant rescattering is a prerequisite for thermalization, 
it is imperative to demonstrate an understanding of the 
centrality dependence of the AA results in order to un- 
derstand the effects of rescattering. While those effects 
may be minor on extensive observables, like the particle 
multiplicity or transverse energy, they should be most 
evident on the momentum distribution of the produced 
particles. Recently it has been argued that a parton cas- 
cade description could successfully describe many of the 
features of central Pb+Pb collisions at SPS energies [Q. 
Surprisingly, low momentum transfer soft parton colli- 
sions were found to have little influence on the final ob- 
servables. Similarly, recent perturbative QCD calcula- 
tions were able to reproduce the preliminary WA98 neu- 
tral pion result for central collisions without need 
for the effects of parton energy loss or rescattering Q . In 
this letter we present neutral pion spectra for 158-A GeV 
208p|^_l_208p|-| collisions and investigate in detail the cen- 
trality dependence of the spectral shape and yield. 
The CERN experiment WA98 ||,| consists of large 



acceptance photon and hadron spectrometers together 
with several other large acceptance devices which allow 
to measure various global variables on an event-by-event 
basis. The results presented here were obtained from 
an analysis of the data taken with Pb beams in 1995 and 
1996. The minimum bias reactions {(JminMas ~ 6300 mb) 
are divided into eight centrality classes using the trans- 
verse energy Et measured in the MIR AC calorimeter. In 
total, K, 9.6 • 10^ reactions have been analyzed. 

Neutral pions are reconstructed via their 77 de- 
cay branch using the WA98 lead-glass photon detector, 
LEDA, which consisted of 10,080 individual modules with 
photomultiplier readout. The detector was located at a 
distance of 21.5 m from the target and covered the pseu- 
dorapidity interval 2.35 < 77 < 2.95. 

The general analysis procedure is similar to that used 
in the WA80 experiment and described in . Hits in the 
lead-glass detector are combined in pairs to provide dis- 
tributions of pair mass vs. pair transverse momentum (or 
transverse mass) for all possible combinations. Subtrac- 
tion of the combinatorial background is performed using 
mixed event distributions. The resulting momentum dis- 
tributions are corrected for geometrical acceptance and 
reconstruction efficiency. The efficiency depends on the 
particle occupancy in the detector and therefore has been 
calculated independently for each centrality bin. The sys- 
tematic error of the pion yields is mainly due to errors 
in the reconstruction efficiency for central collisions and 
to corrections for non-target interactions for peripheral 
collisions. The systematic error on the absolute yield 
is « 10% and increases sharply below pT = 0.4GeV/c. 
An additional systematic error originates from the un- 
certainty of the momentum scale of 1 %. The influence 
of this rises slowly for higher pT and leads to an error of 
15 % at Pt = 4 GeV/c. A detailed discussion of the anal- 
ysis procedure and the error contributions will be given 
in a forthcoming publication. 

The measured neutral pion spectrum from central 
Pb+Pb reactions (10% of min.bias cross section) as a 
function of ttit — tjq is shown in Fig. ^ The data are 
compared to predictions of the string model Monte Carlo 
generators FRITIOF 7.02 ^ and VENUS 4.12 0. As 
already observed in S-|-Au reactions both generators 
fail to describe the data well at large m^. The FRITIOF 
prediction is more than an order of magnitude lower 
at high niT while VENUS significantly overpredicts the 
data. Alternatively, it has recently been shown that per- 
turbative QCD calculations, including initial state mul- 
tiple scattering and intrinsic pT , are able to describe 
the preliminary WA98 data at intermediate and high px- 
This prediction is included in Fig. |l| as a solid lincQ The 
pQCD calculation shows a very good agreement in the 
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high rriT region. This surprising agreement has been in- 
terpreted as an indication for unexpectedly small effects 
of parton energy loss 0. On the other hand, the par- 
ton cascade Monte Carlo code, VNI, which provides a 
more detailed pQCD description, overpredicts the mea- 
sured WA98 result by more than a factor of ten at large 

An alternative picture of the ultra-relativistic heavy 
ion reaction assumes that the mean-free path of the pro- 
duced particles in the overlap region is sufficiently small 
that a thermal or hydrodynamical interpretation may be 
applied. The excited matter is characterized by thermal 
variables, such as temperature and collective velocities, 
which evolve over time. Hydrodynamical descriptions 
(see e.g. [|2|) with properly adjusted parameters can de- 
scribe the momentum spectra reasonably well. This is, 
however, beyond the scope of this paper - a detailed 
investigation of hydrodynamical descriptions of the mo- 
mentum spectra will be given in a forthcoming publica- 
tion. 

In view of the above discussion and the difficulty 
to describe the details of the neutral pion spectrum, 
it is apparent that the theoretical description of ultra- 
relativistic nucleus-nucleus collisions remains uncertain. 
In order to demonstrate a consistent description of nu- 
clear effects it is important to investigate the details of 
the pion production as a function of the system size. To 
study the centrality dependence of the spectral shape in 
a manner which is independent of model or fit function 
we have used the truncated mean transverse momentum 
{Pt{pT"))^ where 

(PTipm) = / PT^dpT / / -^dpT -pT 

\Jp,^^^ dpT I Jp™- dpT ) 

(1) 

The lower cutoff p™"' = 0.4 GeV/c is introduced to avoid 
systematic errors from extrapolation to low and has 
been chosen according to the lowest pt of the present 
data where systematic uncertainties imposed by the nec- 
essary corrections are still small. In general, the value 
of {pt{Pt™)) differs from the true average pr, except 
in the case of a purely exponential distribution da/dpx- 
For a purely exponential invariant cross section, d^a /dp^, 
{pt (p™™ ) ) decreases with increasing p™™ . 

Figure g shows (pt(Pt*")) as a function of the av- 
erage number of participants Npart for 158-A GeV 
208pj^_l_pj^ collisions. For comparison, (pt{Pt^^)) 'val- 
ues for 200-A GeV S+Au Q and from a parametrization 
of pp data are also included. Npart is extracted by 



error by the author of |7| and have changed by ~ 10 — 30% 
compared to the publication. 



the assumption of a monotonic relation between impact 
parameter and transverse energy and using the result- 
ing correspondence between measured cross section and 
impact parameter. The average number of participants 
is calculated from nuclear geometry using the extracted 
impact parameter. Together these data show the gen- 
eral trend of a rapid increase of {pt{p™^)) compared 
to pp results for small system sizes. For Npart greater 
than about 30 the mean transverse momentum appears 
to attain a limiting value of « 280MeV/c. Venus 4.12 
pT| calculations show a qualitatively similar behaviour, 
although the values of (pt(Pt™)) ^^"^ somewhat lower 
than the experimental data. The simple implementation 
of rescattering which is used in this model seems to be 
strong enough to lead to a saturation for semi-peripheral 
collisions as in the experimental data. One should, how- 
ever, keep in mind that Venus 4.12 does not correctly 
describe pion production at high pT (see figure |l|) . 

Alternatively, exponential fits to the spectra over the 
range 0.7GcV/c^ < rax — mo < 1.9GeV/c^ have been 
performed. Although they do not adequately describe 
the spectra, even in this limited transverse mass range, 
they provide another measure of the spectral shape. The 
fitted inverse slopes show a similar saturation effect . 

Earlier investigations of the dependence of (px) of pi- 
ous on system size |^,|l^,^ at SPS energies have sug- 
gested such a saturation for large systems. The present 
study is the first investigation of the dependence with Pb 
ions at the SPS. Preliminary results from the AGS have 
indicated a weak increase in the average rriT of pious with 
the number of participants for Au-t-Au collisions p^ . 

It is important to note that the observed limiting be- 
haviour is very different from the observations in pp or 
pp collisions. For very high energy (px) rises with the 
pseudorapidity density of charged particles jl^-^^. In 
that case, more violent parton scatterings presumably 
result in a harder spectrum of leading particles together 
with a greater multiplicity of fragmentation products. 
This would lead to the observed correlation between (px) 
and multiplicity. At lower ^/s, comparable to the data 
presented here, (px) decreases for increasing multiplicity 
pl| , most likely due to energy conservation. In the case 
of nuclear reactions, this anti-correlation is lost due to the 
large number of binary collisions. Instead, the initial in- 
crease of (pt(Pt™)) with Npart is interpreted as a result 
of multiple scattering. Initial state multiple scattering, as 
suggested as explanation for the Cronin effect [Q , would 
imply a continuing increase of (pt(p™*")) for more cen- 
tral collisions. Here however, the surprising observation 
is that additional multiple scattering, implied by increas- 
ing Npart, does not alter the pion distributions. This is 
most easily understood as a consequence of final state 
rescattering and is, of course, the behaviour expected for 
a thermalized system. 

More detailed information about the centrality depen- 
dence of the pion spectral shape and yield is shown in 
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Fig. ^ where the neutral pion yield per event has been 
parameterized as E(PN/dp^ oc Np^^f^ ■ (To(pt)- The re- 
sults for Npart > 30 are well described by this scaling 
with an exponent aipr) ~ 1-3, independent oi pt- Con- 
sistent with the previous discussion, the results indicate 
a constant spectral shape over the entire interval of mea- 
surement from 0.5 < PT < 3 GeV/c. The observed Npi% 
scaling for symmetric systems implies a scaling with the 
number of nucleon collisions, as confirmed by a simi- 
lar analysis. However, this scaling does not extrapolate 
from the pp results. On the contrary, when comparing 
semi-peripheral Pb-f Pb collisions with pp the exponent 
a varies over the entire pt interval, confirming the very 
different spectral shapes. 

In summary, we have analyzed the centrality depen- 
dence of high precision transverse momentum spectra of 
neutral pions from 158AGeV Pb-|-Pb collisions. The neu- 
tral pion spectra are observed to show increasing devia- 
tion from pp results with increasing centrality, indicating 
the importance of multiple scattering effects. However, 
for centralities with more than about 30 participating nu- 
cleons, the shape of the transverse momentum spectrum 
becomes invariant over the interval 0.5 < pt < 3 GeV/c. 
In this interval the pion yield scales like N^^^^, or like 
the number of nucleon collisions, for this range of cen- 
tralities. Since the amount of rescattcring increases with 
centrality, the invariance of the spectral shape with re- 
spect to the number of rescatterings, most naturally sug- 
gests a dominantly thermal emission process. It will be 
important to determine whether cascade models which 
reproduce the observed invariant spectral shape will sup- 
port the interpretation as an "effective" thermalization 
due to significant rescattcring. 
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FIG. 1. Transverse mass spectra of neutral pions in central 
collisions of 158 AGeV Pb+Pb. Invariant yields per event 
are compared to calculations using the FRITIOF 7.02 
and VENUS 4.12 jnl Monte Carlo programs. Predictions of 
a pQCD calculation|7| are included as a solid line. The inset 
shows the ratios of the results of the Monte Carlo codes to 
the experimental data. 
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FIG. 2. Truncated mean transverse momentum (pt(pt'")) 
of tt" mesons as defined by Eq. ^ plotted as a function of the 
average number of participants Npart ■ The solid circles corre- 
spond to the 8 Et based centrality selections for Pb+Pb. The 
open square shows {pt(Pt^^)) extracted from a parametriza- 
tion of pp data scaled to the same cms-energy llSl, the open 
circles the results for S+Au collisions at 200 ylGeV For 
comparison, results from Venus 4.12 ||ll| are included as his- 
tograms for Pb+Pb collisions and as a star for pp. A cut 
parameter py'" = 0.4GeV/c was used. 
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FIG. 3. The exponent a{pT) of the dependence of the tt" 
yield on the average number of participants Npart plotted as a 
function of the transverse momentum for 158 AGeV Pb+Pb. 
The solid circles are calculated based on the centrality selec- 
tions with Npart > 30. The open circles are calculated based 
on the ratio of the semi-peripheral data (Npart ~ 45) to a 
parameterization of pp data. 
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